Free volume distributions inside a bidimensional granular medium 
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We investigate experimentally within a two-dimensional cell the distribution of the free volume associated either to single 
grains or to clusters of grains. Our main result is that the logarithm of the probability to find a free volume per grain 



N in a cluster of N grains scales as N a g(v J N ), with a = 0.75. We interpret this non extensive scaling factor N a as 
an evidence for the onset of long range correlations between the free volumes of individual grains. We also discuss the 
possible relation between g{v^ N ) and Edwards entropy. 



1 INTRODUCTION 

Everyday life and recent experiments sug- 
gest that a thermodynamical descrip- 
tion of granular media might be feasi- 
ble (fN owak, Knight, Ben-Naim, Jae ger, and Nagel 1998} 
|Schroter, G oldman, and Swinney 2005). Given that gran- 
ular media consist in a large number of grains, there is a 
strong motivation for providing a statistical ground to this 
hypothetic thermodynamical description. However, the 
stationary dynamics of granular media results from the 
balance of dissipation and forcing by a non thermal source, 
so that it does not leave any known ensemble invariant 
and little is kn own on the typical configurations explore d 
dynamically (Barrat, Kurcha n, Loreto, and Sellitto 2 000). 

Still, it has been argued by Edwards and 

collaborators dEdwards and O akeshott 19891 

Meht a and Edwards 1989ft that the dynamics is con- 
troled by the mechanically stable -the so-called blocked- 
configurations and that all such configurations of a 
given volume are statistically equivalent provided that 
the driving involves extensive manipulations, such as 
shaking, shearing or pouring. This immediately leads 
to the definition of a configurational entropy S con f(V), 
and the associated state variable, the "compactivity" 
Klf = 9S conf /dV. 

First attempts to test this flat measure assumption 
have taken advantage of the analogy between the gentle 
comp action of grains and the aging of glassy sys - 
tems ( |Liu and Nagel 1998| ID' Anna and Gremaud2 001 ). 
Despite some cl ear examples where Ed wards' ap- 
proach fails (Godreche and Luck 2005), explicit 
checks have been made so far in mean field mod- 
els of the glass transition (Monasson 1995J), in 
schematic finite-dimensional models with k inetic 
constraints (Barrat, Kurchan, Loreto, and Sellitto 2000), 



in spin glass models with n on-thermal driving be- 
tween the blocked states (Dea n and Lefevre 2001ft 
and finally in a f ew more realistic models of parti- 
cle deposition ( Bre y, Prados, and Sanchez-Rey 2000 ) 
or MP s imulatio ns of she ar driven granular me- 
dia (IMakse and Kurchan 2003 . 

However, clear evidence in real granular media is 
still lacking. A steady state, history independent dy- 
namics un der vibration has been observed by Nowak e t 
al. ( |Nowak, Knight, Ben-Naim, Jaeger, and Na gel 1998), 
for tapped glass beads at rather high volume fraction. 
More recently, a one to one correspondance between the 
volume fraction fluctuations and the volume fraction it- 
self have been observed, so that assuming the valid- 
ity of the usual thermodynamical relation < 5V 2 >= 
X 2 d < V > /dX, a measure of the compactivity X 
as a function of the volume fraction ha s been de- 
rived ( Schroter, Goldman, and Swinney 2005 ). 

It would now be of major interest to produce further ex- 
perimental evidence of the statistical foundation of such 
thermodynamical properties. In the present paper, we re- 
port experimental results on the free volume distributions 
inside a bidimensional granular medium. We show that the 
probability distribution of the free volume per grain inside 
a cluster of N grains follows: 
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2 EXPERIMENTAL SET-UP 

The experimental cell (Fig. [It) consists in a 60cm by 50cm 
rectangular glass container with an inside gap of thick- 
ness 3.5mm. It contains 5000 nickel plated brass cylin- 
drical spacers of two different sizes in equal number. Both 
kind of spacer is 3 mm thick, and their diameters are re- 
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spectively d s = 4mm and d\ = 5mm. In the following d s 
has been chosen as the unit length, so that the respective 
area of the small and large disks are v® = 7r/4 ~ 0.78 and 
vf = 7r/4(di/d s ) 2 ~ 1.23. The cell is half filled with a sin- 
gle layer of such hard disks mixed together, resulting in an 
homogenous and disordered bidimensional packing. The 
cell is mounted on a horizontal axis and rotated around 
this axis in such a way that the grains fall from one side to 
the other every half cycle. 




Figure 1 : Experimental set-up and sketch of the modified 
Voronoi tessellation. 

The experimental procedure is the following. The cell 
starts in a vertical position and is rotated of one cycle, at a 
constant speed of one cycle per minute. During this cycle, 
the grains fall from one side to the other and then back 
to the initial side. We then stop the engine, wait for a few 
seconds to allow the system to reach a mechanically stable 
state, and take a picture of the bulk with a CCD camera. 
We perform 15 000 of such cycles. 

The pictures hence taken display on average 300 grains. 
For each picture the centers of the spacers are located and 
their Voronoi diagram is computed (Fig [II)), taking into 
account the bidispersity of the assembly. We then collect 
the area of the cells along with the type, position and index 
of the associated grains. Out of these raw data, we extract 
and analyze the statistical distribution of the volumes oc- 
cupied first by one grain, then by clusters of an increasing 
number of neighboring grains. 

3 DATA ANALYSIS 
3.1 One grain volume distribution 
Figure displays the distribution of the Voronoi cell area. 
The distribution displays two peaks centered on < v s >= 
1.00 (resp. < v\ >— 1.49), the averaged area occupied 
by the small, (resp. the large) grains computed indepen- 
dently. Also indicated on the figure are the minimal values 
that a Voronoi cell can possibly take - the closest regular 
hexagon- for each type of grain, v™ in = y/3/2 ~ 0.866 
and vp in = V3/2(di/d s ) 2 ~ 1.35. Both peaks present a 
well defined exponential tail, which is easily isolated when 
considering the distributions of the free volume (v^ sl = 
v — v™{ n ), for each type of grain as shown on the inset of 
figure |5| It shall be noticed that the characteristic free vol- 



are closer to each other than one would have expected 
given the particule size ratio. In other words, the character- 
istic free volume accessible to one grain not only depends 
on the size of that grain, but also on the neighborhood, 
most often composed of grains of both sizes. 




Figure 2: Distribution of the Voronoi cells area. Vertical 
dashed lines : minimal Voronoi cell area. Vertical dash- 
dotted lines : conditional average Voronoi cell area. Inset: 
distributions of the free volume conditioned by the grain 
size; (dark): small grains; (grey): large grains 
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Figure 3: Distribution of the free volume including both 
kinds of grains. Inset: density of states of one grain sur- 
rounded by a free volume (see text for details). 

Figure [3 displays the distribution of the free volume 
without conditioning it by the type of grains. One still ob- 
serves the well defined exponential tail, with a character- 
istic free volume = 0.058. In the spirit of Edwards' 
description, one may assimilate this volume, characteris- 
tic of the sample, to a compactivity. This is of course a 
rather strong interpretation of the data, since it would be 
necessary at least to check that the dependence of the dis- 
tributions on this compactivity is indeed fully embedded 
in the exponential tail. Still, as an indication, we have 
plotted in the inset of figure [3] the quantity \n(p(v^)) — 
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\n.(P(vf)) + vf /vf*, the density of states of one grain 
surrounded by a free volume . It saturates for large free 
volumes, after a sharp increase for small values. 

At this point, it should be noticed that in granular me- 
dia, one may not expect the N-body probability distribu- 
tion for the volumes of a set of N grains to factorize as 
a product of one-grain distributions, because of possible 
correlations and degeneracies of the free volume states. In 
order to check if one recovers the neat properties of the 
usual canonical description for a large enough subsystem 
of the grain assembly, we hence turn to the volume distri- 
butions for clusters of neighboring grains. 

3.2 Free volume distribution for clusters of neighboring 
grains 

Now dealing with clusters of grains, the free volume and 
the full volume are no longer equivalent, since clusters are 
composed of both kinds of grain. Figure |4] (respectively 
figure 0> displays the distribution of the volume (resp. the 
free volume) per grain inside clusters of N neighboring 
grains, v N = iV^E^-^ (resp. v N 7V" 1 H^ 1 v/ ). 
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Figure 4: Distribution of the volume per grain inside clus- 
ters of N grains: the larger N, the narrower the distribution 

The distribution of the free volume per grain v N con- 
verges much faster towards a well defined distribution than 
that of the volume. Indeed, it is clear from figure 0J in the 
case of the volume distributions, as N increases from 1 
to 10, that combinatory effects induced by the bidispersity 
give rise to series of peaks which complexify the distri- 
bution. For N = 2, one observes three peaks, which are 
related to the three two-particles cluster configurations. 
Similarly, for TV = 3 -not shown here-, the distribution 
exhibits four peaks, etc. Such effects are absent in the case 
of the free volume distributions (figure |5j. Also, describ- 
ing the volume distribution imposes to know the minimal 
value of the volume per grain, a non trivial quantity for 
clusters composed of grains of both sizes, whereas by con- 
struction, the distribution of the free volume readily ad- 
dresses this issue. Accordingly, we restrict the following 
analysis to the free volume distributions. 

We propose to describe the distributions of the free vol- 
ume per grain inside a cluster of N neighboring grains by 




Figure 5: Distribution of the free volume per grain inside 
clusters of N grains: the larger N, the narrower the distri- 
bution 



a Gamma law of parameters tjn and Xn : 
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where T is the Euler Gamma function. This choice is 
motivated by the following considerations. In the hypo- 
thetic case where the one-grain free volume distribution 
would follow a Gamma law itself of parameter (rji,Xi), 
if the free volumes inside a cluster were independent 
variables, one would have exactly the proposed Gamma 
law for the distributions of v N , together with the rela- 
tions (t]n = Ntji,Xn — Xi/N). This is a standard prop- 
erty, which seemingly extends for large enough N in the 
case where the one-grain distributio n is not a Gamma 
law dBertin, Dauchot, and Droz2 QQ4). Given that in the 
present case, the general shape of the one-grain free vol- 
ume distribution already shares some of the key property 
of the gamma-law -being defined for > and exhibit- 
ing an exponential dependence for large v$-\ given that 
the distributions seemingly belong to the same functional 
family irrespectively of N, it is reasonable to believe that 
the convergence towards the gaussian law expected in the 
limit of very large N will occur via Gamma laws. 

Once chosen the form of the distributions, we compute 
their first two moments and obtain t]n and Xn, through 



the relations < v J N > 
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f 

t]nX n and < v J N > 
J 



< v f N > 2 = 



t]nXn • As expected, < v J N > rapidly evolves towards 



> 



< 



a constant (figure 0-left). On the contrary < v 

v f N > 2 varies like N~ a with a = 0.75 ± 0.0025, in con- 
trast with the 1/N dependence expected for independent 
variables (figure fright). 

As a direct consequence of the above analysis, the dis- 
tribution of the free volume per grain inside clusters of 
N grains is well described by a Gamma law, the param- 
eters of which exhibit the following dependences on N : 
f] N = VeffN a and Xn = X ef fN~ a , with r] = 3.5 and 
X = 0.041. We have also plotted on these figures the de- 
pendences extracted from direct fit of the distributions by 
the proposed gamma-law; the results are the same. 
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Figure 6: Dependence on N of the first (left) and second 
(right) moments of the free volume distributions (o): com- 
puted from the data; (*): extracted from the fit of the dis- 
tributions by a Gamma law; (plain line): fit of their depen- 
dence on N 

4 DISCUSSION AND CONCLUSION 
Altogether, rewriting the above Gamma law in the limit 
of large N, one obtains that the logarithm of the dis- 
tribution of the free volume per grain inside clusters 
of N grains scales as N a g(v,r] e ff J X e ff) with g(v) — 
r](ln(v/(r]X)) - v/(j)X) + 1), a ~ 3/4, rj eff ~ 7/2 and 
X e ff = 0.041. This central result deserve a few com- 
ments. 

In the case of Poisson Voronoi diagrams in two 
dimensions -where the centers of the Voronoi cells 
are rando mly chosen-, it has be en shown very re- 
cently (IJarai-Szabo and Neda 2004 ) that the distribution 
of the Voronoi cell areas normalized by the mean area 
is precisely well described by a Gamma law of parame- 
ter r] = 7/2 and X = 1/rj. In the absence of correlations 
one obtains in the limit of large N, that the logarithm of 
the distribution of the volume per cell inside clusters of TV 
grains scales as Ng(v, 77, X). Let us now compare our ex- 
perimental results to this academic situations of ponctual 
grains with no hard sphere exclusion, and no correlations. 

First, the observed non extensive factor N a is presum- 
ably the evidence of long range correlations between the 
free volumes of individual grains. Indeed, in the presence 
of correlations decaying with the distance r as 1/r 7 , one 
has in two dimensions, for 7 < 2, that the second moment 
of the average of N centered random variable scales like 
N~ 7 / 2 . In the present case, we would thus infer the exis- 
tence of long range correlations decaying like 1 / r 3 / 2 . 

Second the above analysis has allowed us to define ef- 
fective parameters for the probability distribution of the 
free volume for one grain. Rather surprisingly rj e ff = 77: 
it seems that the hard sphere exclusion for the grains does 
not constraint the value of this parameter. 

Finally, one can also write the distribution of the free 
volume per grain inside clusters of N grains as: 

P(v) = ^-e- Na (^- siv) ), with s (v)=r](ln^- + l), 

Z T]X 



and thereby x = §f 



(v) 



, an exact result given the Gamma 



law distribution and more generally expected from a sad- 
dle point calculation in the large N limit. 

Future work will have to check the robustness of the 
above results, when varying the type of grains or the dy- 
namical protocole. Also, it will be of major interest to be 



able to vary the averaged volume fraction of the sample, in 
order to investigate the equation of state in the light of the 
above statistical description. 
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